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THE  RESPONSE  OF  MECHANICAL  SYSTEMS  TO  RANDOM  EXCITATION 


It  ^  oduction 


Mechanical  aystems  are  not  always  excited  by  a  harmonic  force  of 
fixed  frequency  and  amplit\ide.  Often  the  excitation  input  is  of  random 
nature,  and  the  response  of  the  system  displays  no  orderly  trends.  Instan¬ 
taneous  values  and  phase  are  meaningless  in  such  cases,  and  the  problem 
must  be  treated  from  a  statistical  approach.  It  is  the  purpose  of  this 
paper  to  outline  such  an  approach  as  related  to  the  dynamic  response  of 
structures. 

The  mathematical  basis  of  the  statistical  technique  has  been  exhaus¬ 
tively  treated  by  Rice.^  Engineering  applications  of  statistical  concepts 
have  been  followed  with  notable  advances  in  fields  such  as  communication. 
Applications  to  the  structural  field  have  only  recently  received  attention, 
with  treatment  of  problems  of  naturally  statistical  natxire  such  as  buffeting*^ 
and  fatigue.'^  With  greater  utiliiation  of  rocket  and  Jet  engines,  the  statis¬ 
tical  approach  to  structural  dynamics  asstxmes  a  role  of  increasing  impor¬ 
tance. 

Two  problems  of  interest  to  missile  design  have  been  treated  here. 

The  first  problem  deals  with  the  longitudinal  vibration  of  a  slender  rod 
excited  by  a  random  force  at  one  end.  The  second  problem  is  that  of 
flexural  vibration  of  a  free-free  beam  excited  by  a  random  transverse 
force  at  one  end.  Briefly,  the  problem  require#  statistical  description 
of  the  random  excitation  and  determination  of  the  frequency  reSf>onse  of 
the  structure  including  structural  damping.  With  this  information,  it 
IS  possible  to  determine  i'or  any  point  in  the  structure  the  probability  of 
exceeding  any  specified  response  such  as  stress,  deflection,  moment, 
or  other  quantities  of  interest. 
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Fundamental  Concepts 


I 


To  establish  certain  fundamental  concepts,  we  will  consider  a  linear 
system  of  single  degree  of  freedom  excited  first  by  a  harmonic  force 
P  cos  ut.  The  differential  equation  for  such  a  system  may  be  written  as 


y  +  Zr  y  + 
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(1) 


where  only  the  real  part  of  this  equation  is  considered.  Its  steady- state 
solution  is  the  real  part  of  the  equation. 
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which  upon  substitution  into  equation  (1)  leads  to  the  well-known  resvilts 
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The  quantity  of  equation  (3)  is  the  impedance  function  defined  as  the 
ratio  of  the  input  to  the  output.  Impedance  in  a  mechanical  system  is 
usually  defined  as  the  input  force  divided  by  the  output  velocity;  however, 
for  this  discussion,  we  will  apply  the  concept  in  a  more  general  sense  of 
input  divtfed  by  the  output  without  designating  the  quantities  involved. 

Of  interest  here  is  the  mean  square-response  defined  by  the  equation 
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If  now  the  frequencies  ujj,  '  •  <^iffer  only  by  a  small  amount,  we 

approach  a  continuous  spectrum,  and  the  mean  square-response  in  the 
frequency  interval  Aw  becomes 
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Hence,  by  comparison  with  eqtiation  (9),  we  arrive  at  the  resvilt 
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Equation  (14)  is  the  general  equation  for  the  mean-square  response  of 
the  structure  in  terms  of  the  power  spectral  density  of  the  input  and 
the  impedance  function  of  the  system. 

Evaluation  of  y^  for  Systems  with  Small  Damping 


If  the  damping,'^  ,  is  small,  Z(w)  will  undergo  a  large  change  near 
the  resonant  frequency,  w^;  and,  if  the  variation  in  f(w)  is  of  lesser 
extent  in  this  neighborhood,  equation  (14)  can  be  approximated  with 
good  accuracy  by  the  expression 
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To  gain  some  useful  concepts  regarding  this  expression,  we  will  consider 
the  single-degree-of-freedom  system  of  the  previous  section  where, 
omitting  the  factor  k^,  the  impedance  equation  is 
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If  now  the  frequencie*  Qj,  only  by  a  small  amount,  we 

approach  a  continuous  spectrum,  and  the  mean  square-response  in  the 
frequency  interval  Aw  becomes 
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Equation  (14)  is  the  general  eqviation  for  the  mean-square  response  of 
the  structure  in  terms  of  the  power  spectral  density  of  the  input  and 
the  impedance  function  of  the  system. 

Evaluation  of  y^  for  Systems  with  Small  Damping 


If  the  damping,^  ,  is  small,  Z(w)  will  undergo  a  large  change  near 

t 

the  resonant  frequency,  w^;  and,  if  the  variation  in  f(w)  is  of  lesser 
extent  in  this  neighborhood,  equation  (14)  can  be  approximated  with 
good  accuracy  by  the  expression 
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To  gain  some  useful  concepts  regarding  this  expression,  we  will  consider 
the  single-degree-of-freedom  system  of  the  previous  section  where, 
omitting  the  factor  k2,  the  impedance  equation  is 
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A  plot  of  eqxi&tion  ( i6),  shown  in  Figure  1,  indicates  its  peak  value  to 
be  1/4^^  at  =  1-  It  can  be  shown  by  solving  the  equatioa 
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for  small  that  the  frequencies  at  the  half-power  level  are 


which  determines  the  width  of  the  lA|Z(w){  curve  at  thia  point  to  be 

..ft:  .  .  2 

The  area  iinder  the  1/ jZ(w)|  curve  for  small  ^  can  be  evaluated 
'  by  letting  =  1  +  c.  where  c  <<  1.  The  integral  of  equation  ( 1 5) 
'i  for  small  t  then  becomes 
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Since  t  is  assiimed  small,  c  can  be  several  times  Z  and  still  remain 
small  compared  to  unity,  and  hence  this  integral  is  given  with  sufficient 
accuracy  by 


I  =  ?T  "n- 

The  mean-square  response  then  becomes 

=  H-a»,  f(u  ).  (21 

It  sho\ild  be  noted  here  that  the  integrated  area  is  equal  to  t/2  times  the 
peak  val\je  times  the  width  at  the  half-power  point,  a  relationship  which 
vkull  be  of  use  later. 
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Continuous  and  Multidegree  Freedom  Systems 


Although  equation  (21)  was  derived  from  Z(u))  associated  with  a 
system  of  single  degree  of  freedom,  it  can  be  shown  to  be  applicable 
to  continuous  and  multidegree  freedom  systems  with  slight  modification. 
In  such  cases,  l/|Z(u))|  has  many  peaks;  however,  the  integral 

du,/|z(w)|  at  each  peak  has  the  same  form  as  equation  (19)  except 
for  a  term  which  is  a  function  of  position  x  in  the  system  and  which  is 
not  involved  in  the  evalviation  of  the  integral. 

The  impedance  function  in  the  neighborhood  of  resonance  can  here 
be  expressed  in  the  form 
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where  K{x,  )  I r  is  the  peak  value  and  r  is  the  term  associated  with 

damping.  Again  the  integral  du>/iZ(cj)j  can  be  eval\iated  by  the  rule 
of  the  previous  section,  and  the  mean-square  response  for  such  systems 
take  the  form 


K(x.  w  )  f(w  ).  (23) 
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Illustrative  examples  in  the  following  sections  will  clarify  this  equation. 


Longitudinal  Motion  of  Slender  Rods  under  Random  Excitation 

We  will  consider  here  a  slender  rod  excited  axially  at  the  end, 
x  =  0,  by  a  random  force,  F(t),  and  free  at  the  other  end,  x  =  £. 
Structural  damping  will  be  included  by  a  complex  stiffness,  E(1  +  ig). 
The  differential  equation  of  motion  for  an  element  of  the  bar  is 
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where 
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c  =  y  ~  velocity  of  the  compresaional  w»ve 


=  percentage  of  structural  damping  (in  decimals). 


Taking  the  Laplace  transformation  and  making  the  subsitution 


(1 1  ig)"'  =  (I  -ig). 


the  subsidiary  equation  becomes  > 
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The  general  sclution  for  the  above  equation  is 


u  (x,  a)  =  C, 
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Fitting  this  solution  to  the  boundary  conditions. 
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we  arrive  at  an  expression  for  the  subsidiary  stress,  s  (x,  s),  which  is 
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To  determine  1/  |Z(w)j^,  we  evaluate  this  equation  for  a  harmonic 
excitation,  F(o,  t)  =  F^  which  results  in  the  expression 
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If  we  define  the  impedance,  jz(b))|,  for  this  problem  to  be  the 
ratio  of  the  input  stress  to  the  output  stress  at  any  position  x,  then 
l/|Z(u)j  is  given  by  the  absolute  value  of  the  quantity  in  the  braces 
in  equation  (30).  Making  the  substitution  \/l  —  ig  =  1  —  i(g/2)  where 
g  is  small,  we  obtain  the  desired  result  to  be 


We  note  in  this  expression  that  if  the  dampijog  g  =  0.  tJien  l/|z(w)|^ 
becomes  infinite  when 
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thus  identifying  the  natural  frequencies  of  the  bar  to  be 
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Retaining  g,  the  impedance  ftinction  in  the  neighborhood  of  resonance 
becomes 
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which  is  in  the  form  of  equation  (22),  where  r 
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Thus  by  substitution  into  equation  (23),  or  by  using  the  previous  area 
rule,  the  mean- square  response  becomes 
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which  can  be  evalviated  when  the  power  spectral  density,  f(u),  of  the 
input  stress  is  known. 


Power  Spectral  Density  of  Excitation 

The  curve  for  the  power  spectral  density,  f(fai),  may  have  variotis 
shapes  depending  on  the  nature  of  the  source  of  excitation.  However, 
from  energy  considerations,  f(u)  for  all  forma  of  excitation  must 
approach  zero  as  u  approaches  infinity.  In  all  cases  f((,»)  maybe 
defined  by  its  mean-square  value  and  some  characteristic  frequency, 


(34) 


One  variation  of  f  (u)  frequently  considered  is  the  monotonically 
decreasing  curve  of  Figure  2,  defined  by  the  equation 
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It  should  be  noted  that  the  above  equation  aatisfie  the  original  ' 
definition 
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The  choice  of  the  characteristic  frequency,  w^,  is  arbitrary; 
however,  it  is  convenient  to  relate  it  to  the  frequency  corresponding 
to  the  half- power  point  of  the  input  spectral  density.  If  repre¬ 
sents  this  half-power  point,  then  u)^  =  2.  1  84%  of  the  input 

energy  is  contained  in  the  frequency  range  0  to 


Another  variation  of  f(b>)  which  is  of  interest  is  the  rectangular 
distribution  of  Figure  3  with  a  cut-off  frequency,  The  spectral 

density,  f(w).  is  then  defined  by  the  equation 
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This  distribution,  v/hich  is  often  referred  to  as  "white  noise,"  can  some¬ 
times  be  used  as  an  equivalent  spectrum  for  the  more  general  types  of 
distribution. 


To  establish  certain  concepts  regarding  the  excitation  spectrom.  it 
is  well  to  discuss  how  f(u»)  is  determined  by  experiment-  The  essential 
components  of  the  measurement  apparatus,  shown  in  Figure  4,  consists 
of  a  band- pass  filter  and  a  meter  or  indicator  which  will  read  the  mean- 
square  values.  Equipment  is  commercially  available  where  the  central 
frequency  of  the  band-pass  filter  is  continually  sw'jpt  and  the  spectrum 
for  f(u)  is  displayed  on  an  oscilloscope.  With  a  continuous  spectrum,  the 
band  of  frequencies  passed  by  the  filter  represents  a  sum  of  the  harmonic 
oscillations  of  frequencies  differing  by  increments,  bu.  Noting  that 
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it  is  evident  that  the  result  is  a  large  number  of  amplitude -modulated 
oscillations  approaching  a  random  amplitude  fluctuation  at  frequencies 
of  order  6(j. 


In  interpreting  equations  (35)  or  (36),  one  must  not  assume  that  the 

height  of  the  f(w)  curve  diminishes  with  large  values  of  the  cut-off  fre- 

quency,  The  spectral  analyzer  in  indicating  f(w)  =  AF  / Au  is 

unaware  of  the  excitation  outside  the  instantaneous  pass  band  A(»»,  and 

hence  f(jj)  is  unaffected  by  the  extent  of  the  frequency  range  of  the 

spectrum.  Thus  the  proper  interpretation  of  these  equations  is  that 

F  / remains  essentially  constant  or  that  the  mean-square  vaJue  of 

the  excitation  increases  with  . 
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Numerical  Evaluation  of  g 

~2 

To  examine  how  the  mean-8q\iare  response,  v  ,  for  the  longitudinal 

motion  of  the  rod  is  related  to  the  mean-square  value  of  the  excitation  stress, 

— 2 

,  we  will  consider  the  two  ty{>e8  of  powder  spectral  densities  discussed 
in  the  previous  section. 

Substitxding  eqxiation  (35)  into  equation  (34),  the  ratio  of  the  mean- 
square  values  for  the  monotonic  spectral  density  is  expressed  by  the 
equation 


In  a  similar  manner,  the  substitution  of  equation  (36)  for  the  rec%- 

tangular  spectrum,  with  an  upoer  limit  on  the  summation  corresponding 

to  the  cut-off  frequency,  w  ,  or  n  =  liw{c/uil),  leads  to  the  result 
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It  is  evident  from  these  equations  that  v  In ^  can  be  plotted  as  a  function 
of  the  nondimensional  quantity,  (c/w^.),  with  x/i,  and  g  as  parameters. 


Results  of  calculations  carried  out  for  x/i  =  1/2  and  g  =  0.01  are 
plotted  in  Figure  5.  These  curves  indicate  that  for  small  values  of 

the  rectangular  and  the  monotonic  spectral  densities  rusult  in 
nearly  the  same  values  of  the  mean-square  ratios,  which  is  not  sur¬ 
prising  when  one  compares  the  above  two  equations.  Since  g  appears 
as  a  linear  factor  in  these  equations,  results  for  other  values  of  damping 
are  obtainable  from  these  curves  by  simple  division. 


Probability  of  Exceeding  a  Specified  Response 


For  a  random  function  such  as  a.  the  normal  probability  distribution 
is  a  reasonable  assumption.  Letting  X  be  any  random  variable  ia  question. 


the  normal  distribution  curve  shown  in  Figure  6  establishes  the  pro¬ 
bability  that  \  will  be  found  in  the  region  X.  to  \  +  dX  to  be 
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(40) 


CX  interest  here  is  the  question  of  the  probability  of  v  exceeding 

some  specified  value,  jt- ,  which  is  n  times  the  root-mean- square 
■—2  ‘ 

responsei^^c  .  To  arrive  at  an  answer  to  this  question,  we  let 
X  =  crl^ZfT^,  in  which  case  the  probability  that  a  lies  between  a  and 
<r  +  dir  is 

dir  (41) 


Thus  the  probability  of  cr  exceeding  a  certain  value,  o-j  = 
is  determined  by  integration  to  be 


(42) 


As  an  example,  if  c/w  i=  0.08,  the  mean-square  ratio  for  the  white- 

~l  — 2 

noise  distribution  is  found  from  Figure  5  to  be  a  / g  =  23.  Thus  the 

*~2  — 2  ^ 

root-mean-square  ratio  is  v  g  =4.80  *  g^, .  The  probability  of  the 
stress  exceeding  some  number,  such  as  3  times  the  root-mean- square 
stress  at  x/j?=  l/2andg  =  0  01,  or  3  time  s  4 .  80  V  g  ^  =  14.4  times 

the  root-mean-square  value  of  the  input  stress,  is 
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Lateral  Motion  of  Beams  Under  Random  Excitation 


■f 

^  We  consider  next  the  lateral  response  of  a  beam  excited  at  the  end, 

^  X  =  0,  by  a  random  force.  F(t),  with  the  end,  x  =  /  ,  free.  Again  we 
I  include  structural  damping  by  a  complex  stiffness,  E(1  +  ig),  and  write 
j  the  differential  eqviation  for  the  beam  element  for  harmonic  excitation. 

^  -TTtW  V  =  0.  ,431 
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m« 

El 

(1  +  ig)’  ^  =  1  -  ig  (44) 

=  ( I  -  ig) 

the  above  equation  becomes 


Letting 


P 


^  =  0. 


(45) 


We  thus  obtain  the  same  solution  as  that  of  the  undamped  beam,  except 
that  p  is  now  replaced  by  its  complex  counterpart,  p.  We  will,  however, 
carry  out  the  solution  in  terms  of  p  and  account  for  p  at  the  very  end. 


Using  the  Laplace  transformation  with  x  as  the  original  variable. 

It  IS  possible  to  arrive  at  a  general  solution  in  terms  of  quantities  at 
X  =  0;  and,  by  successive  differentiation,  equations  for  deflection, 
slope,  moment  and  shear  can  be  obtained  and  expressed  in  matrix  form. 


4 


^  y<*) 

a 

b 

c 

<* 

1  y(0) 

1 

^y(x, 

d 

a 

b 

c 

?  rd)) 

-y  y'Mx) 

c 

d 

a 

b 

-!Ty"(0) 

1 

P 

'  -X  y"’(x) 

b 

c 

d 

a 

-y  y'^MO) 

(46) 


13 


where 


a  =  2*  (cosh  Px  -f  cos  ^x) 

b  =  ^  (sinh  ^x  -f  sin  px) 

c  =  2  (cosh  px  —  cos  px) 


(47) 


d  =  j  (sinh  px  —  sin  px). 

When  X  =  j2,  we  have  the  boundary  conditions  y"(/.)  =  y'"(i)  =  0- 
At  the  end,  x  =  0,  the  boundary  conditions  are  y"(0)  =  0,  and 
l/p^[y'"(0)|  =  P^/P^  El,  where  the  excitation  is  cos  ujt.  With  these 
conditions,  the  remaining  constants  become 
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It  18  now  possible  to  write  the  equation  for  any  quantity  y(x)  to 
y"’(x).  For  instance,  the  equation  for  y"(x)  becomes 
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It  is  evident  here  that  the  natural  frequencies  are  determined  from 
the  equation  D  =  0,  in  which  case  y"(x)  — ►  oo. 


(49) 


Wben  damping  is  included,  we  replace  p  in  the  above  solution  by  p. 
However,  since  small  damping  is  assumed  and  we  are  interested  in  the 
impiedance  function  only  in  the  neighborhood  of  resonance,  we  need  to 
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i 


f 


j 


replace  P  by  p  ozily  in  termi  which  tend  to  sero«  or  in  the  expreeeion 
for  D. 

Letting  p  =*  P  (i—ig)^^^  ^  P(l— i^,  the  expreition  for  D 
becomes 

D  =  (cosh  pi  cos  pi.  —  1)  +  i  ^Pi  (cosh  pi  sin  pi  —  sinh  Pi  cos  pi)  (BO) 

where  terms  containing  g^  have  been  omitted  as  negligible.  We  next 
replace  pi  by  p^i  (1  +  «)  in  the  neighborhood  of  resonance;  and  again 
throwing  ont  iirfimheamTia  ad  hsegLer  es<er  mm  wefil  a* 

cosh  P  I  COM  &  Z~  I  - 
'^n  *^n 

we  arrive  at  the  result 


D  =  (sinh  p^i  cos  p^i  -  cosh  p^  sin  p^i)  (p^i.)  («  -  i  J  ) 


1 

I 


b 

a 


(P^i)(.-i  J). 


(51) 


Substituting  in  this  value  of  D  and  noting  that  in  the  neighborhood  of 
resonance  b(x)  is  negligible  compared  to  the  other  two  terms,  the 
final  equation  for  the  moment  bcconres 


M(x)  =  EIy"(x)  = 


2P  Jl 
o 

T 


{*-  if) 


i  c(x)  + 


c  b 

b  a 


b  d 
a  c 


d(x) 


(52) 


If  we  now  define  the  impedance  as  the  ratio  of  the  input  moment, 

output  moment,  M(x),  the  quantity  of  interest  becomes 


4  -^c(x)  + 


c  b 
b  a 


yB— li 

|a  cl 


d(x)  I 


iz{u)r 


(Pn^) 


(53) 
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By  the  previous  rule  for  integrated  area,  the  mean- square  moment 
contribution  ior  mode  k  is  given  as 


M. 


ir 

7 


^  (l)  ‘ic 

- 

_ 

- 

^(“ll  ,  =  0 

Thus,  summing  over  all  modes,  we  have 


M 


c  b 
b  a 


b  d 
a  c 


Since 


"n  =  (Pn^> 


2  El 


mJ 


and 


(2n+  1)  5 

this  equation  can  also  be  written  as 


772  _  64  \  \l  mi^ 

M  =  — —  -  '  - - -  ■  ■ — 


IT  (2n  +  1) 


c(x)  + 


d(x)  i 


P£  = 


b  d' 
a  c  I 


d(x) 


i. 


M  = 


(54) 


The  squared  quantity  within  the  braces  can  be  identified  as  being  equal 
to  half  the  tabulated  results,  1/2  0  (x),  of  Young  and  Felgar,^  and  for 
this  jroblem  corresponds  to  the  case  of  the  free-free  beam. 


(55) 


(56) 


Numerical  Evaluation  of  M 


As  in  the  longitudinal  case,  we  evaluate  equation  (56)  for  the  two 
types  of  spectral  densities  defined  by  equations  (35)  and  (36).  For  the 
monotonic  variation  of  f(u)),  the  result  reduces  to 
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M  _  128  1  I  El 

—  *  — VT  IT  , 

gw  c  1  ml  ns  1 


U  2 
o 


.  IB 


(2n  +  1)' 


c  b 

c(si  + 

.. 

b  a 

b  d 

a  c 

d(x)  ► 

j 

For  the  white-noise  spectral  density.-.  £(<4  =  /u^,  we  obtain 

the  ratio  of  the  mean  squares  to  be 


n 


M 


M  J_  a  ^  1 

^  **c  \  ml^  n=l  (2n  ^  !)^  1 


M  2 
o 


' 

c 

b 

> 

c(x)  + 

b 

a 

a(xj 

b 

d 

1 

1 

i 

1 

a 

c 

> 

The  summation  in  this  case  is  terminated  at  tlite 

c 

corresponding  to  the  cut-off  frequency,  this  relationship  being 


til  s  (2n  + 
c  c 


El 

ml^ 


<58) 


(59) 


Equations  (57)  and  (58)  were  numerically  evaluated  for  xM  -  i/2 
and  g  s  0.01  and  plotted  in  Figure  7.  These  curves  are  to  be  inter¬ 
preted  in  a  manner  similar  to  that  of  the  longitudinal  case ,  and  the 
probability  equation  (42)  is  again  applicable  with  n  interpreted  as  the 
number  of  times  the  specified  moment  exceeds  the  root -moan- square 
moment. 
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